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ABSTRACT 

The  steady-state  excitation  of  longitudinal  waves  by  a  pair  of 
idealized  grids  immersed  in  a  collisi onlcss  plasma  and  driven  at  a 
frequency  small  compared  with  the  ion  plasma  frequency  is  investigated 
theoretically.  In  linear  theory  the  Pourier-inversion  integral  which 
determines  the  spatial  behavior  of  the  potential  in  the  plasma  is  expressed 
as  a  sum  of  two  integrals  which  embody  the  interactions  of  phase-velocity 
components  of  the  wave  with  ions  and  electrons.  An  appropriate  choice 
of  the  deformed  contour  of  integration  permits  evaluation  of  the  response 
as  the  sum  of  the  residue  of  the  dominant  "ion-acoustic"  pole  and  of  the  two 


branch-cut  integrals,  A  perturbation-series  expansion  of  the  potential  and 
the  species  distribution  functions  in  the  (nonlinear)  Vlasov  equation  yields 
a  hierarchy  of  equations.  In  each  order  the  equations  are  linear  in  the 
perturbation  quantities  of  that  order  and  have  driving  terms  composed  of 
quadratic  combinations  of  lower-order  quantities.  For  sufficiently  small 
amplitude  of  excitation  the  principal  contributions  to  the  response  come 
from  the  first-order  (linearized  Vlasov)  equation  and  the  second-order  equa¬ 
tion.  .  In  second  order  the  steady-state  response  consists  of  zero -frequency 
and  double- frequency  components.  The  second-order  equations  are  Laplace- 
Fourier  transformed  and  resulting  velocity  integrals  are  expressed  in  terms 
o*  plasma  dispersion  functions.  By  approximating  the  driving  terms  by 
their  dominant-pole  component,  one  can  express  the  steady-state  double¬ 
frequency  response  as  a  single  Fouricr-inversion  integral.  As  in  the  linear 
problem,  the  integral  can  be  evaluated  as  the  sum  of  a  residue  and  of  "ion 
like"  and  "elcctron-like"  branch-cut  integrals.  Numerical  results  are 
presented  for  the  linear  and  nonlinear  cases. 
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I.  INTRODUCTION 


The  calculation  by  Gould1  of  the  steady-state  response  in  linear 

theory  of  a  collisionless  plasma  to  grid  excitation  of  longitudinal  waves 

at  frequencies  below  the  ion  plasma  frequency  exhibits  good  agreement 

2 

with  the  experiment  of  Wong,  D’Angelo,  and  Motley,  and  supports  the  inter¬ 
pretation  that  spatial  damping  of  the  wave  is  due  to  Landau  damping.  The 
importance  of  understanding  the  nature  of  the  damping  and  of  studying  non¬ 
linear  processes  in  plasmas  suggest  the  desirability  of  extending  the 
theory  into  the  nonlinear  regime.  As  the  amplitude  of  excitation  is  increased, 
one  expects  the  steady-state  response  to  include  harmonics  of  the  applied 
frequency.  If  the  amplitude  is  sufficiently  small,  only  zero- frequency  and 
double -frequency  harmonics  are  significantly  excited,  and  a  perturbation 
procedure  should  yield  this  nonlinear  response. 

In  Section  II  the  response  in  linear  theory  is  determined  by  a  method 
which  possesses  certain  advantages  over  that  of  Gould.  The  response  is  ex¬ 
pressed  as  the  sum  of  two  branch-cut  integrals  which  contain  the  velocity 
derivatives  of  the  ion  and  the  electron  distribution  functions  and  give  in¬ 
sight  into  the  role  of  the  two  species  in  the  Landau  damping  of  the  wave. 

The  response  may  also  be  expressed  as  the  sum  of  the  residue  of  the  ion- 
acoustic  pole  and  of  the  two  integrals  evaluated  along  deformed  contours 
which  are  well-suited  for  obtaining  high  accuracy.  Calculations  are  per¬ 
formed  both  for  the  case  of  negligible  grid  spacing  (dipole  limit)  considered 
by  Gould  and  for  the  case  of  finite  spacing  between  the  grids;  the  compari¬ 
son  illuminates  the  role  of  the  model  of  the  grid  in  determining  the  char¬ 
acter  of  the  response.  In  Section  III  perturbation* series  expansions  of 
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electric  field  and  species  distribution  functions  are  introduced  into  the 
Vlasov  equation;  the  equations  for  lowcst-order  nonlinear  response  arc  ob¬ 
tained.  Performing  the  appropriate  Laplace-Pourior  transforms,  analyt- 
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ically  continuing  (in  the  manner  of  Landau  )  functions  of  two  complex  vari¬ 
ables  defined  by  velocity  integrals  which  occur  in  the  formulation,  and  re¬ 
lating  these  functions  to  the  plasma  dispersion  function,  a  double  Fourier- 
inversion  integral  for  the  lowcst-order  nonlinear  response  in  the  plasma  is 
obtained.  This  integral  contains  a  quadratic  combination  of  the  electric 
field  in  linear  theory  as  a  driving  term.  In  Section  IV  reduction  of  the 
double  integral  for  the  double-frenuency  component  to  a  single  Fourier-inver- 
sion  integral  is  achieved  by  approximating  the  electric  field  in  linear  theory 
by  its  dominant-pole  component.  The  nonlinear  response  is  expressed  as  the 
sum  of  a  residue  and  of  ion-like  and  electron-like  integrals,  as  in  the 
linear  problem.  In  Section  V  the  branch-cut  integrals  and  the  residue  arc 
considerably  simplified  by  utilizing  the  square  root  of  the  mass  ratio  as  a 
smallness  parameter.  Calculations  are  performed  for  a  range  of  values  of 
the  electron-to-ion  temperature  ratio.  In  Section  VI  the  zero-frequency 
component  of  the  nonlinear  response  is  shown  to  consist  of  a  polarization  of 
the  plasma,  with  no  species  current  densities. 

II.  LINEAR  PROBLEM 

The  steady-state  oscillatory  potential  produced  in  a  uniform  and  in¬ 
finite  collisioniu's  plasma  by  excitation  at  a  frequency  of  a  pair  of 

closely-spaced,  idealized  grids,  which  produces  an  external  oscillating 
charge  density  but  intercepts  no  particles,  is  given  by  <f»(x,t)  =  $(x)exp(-j c^t) 
+  complex  conjugate,  in  which 


Here  is  the  amplitude  of  the  sur face -charge  density  on  either  grid 

and  Xq  is  the  separation  between  the  grids;  the  dipole  limit,  ■+  0, 

0^  x^  =  const,  is  considered,  phe  effect  of  finite  spacing  between  the 
grids  is  examined  below.)  The  "plus"  and  "minus"  dielectric  functions  are 
given  by 


in  which  the  corresponding  plasma  dispersion  functions  are  defined  by 


Z(?> 


-  DO 

exp  dr 

71“  (t  -  5) 


when  Im{^)  is  positive  and  negative,  respectively,  and  by  the  analytic 
continuation  of  these  integrals  elsewhere.  The  functions  Z+(c)  have  a 
branch-cut  along  the  real  axis;  the  situation  may  be  understood  by  consider¬ 
ing  the  integral  to  be  the  limit  of  an  integral  over  the  finite  range  -c  to 
c,  as  c  ■+  «>.  Branch  points  occur  at  t,  =  ±c.  The  mapping  of  the  branch- 
cut  and  branch  points  onto  the  k  plane  for  w  =  Wq  +  ie  is  shown  in  Fig.  la, 
along  with  the  primitive  Fourier  inversion  contour.  The  first  step  in  evalu¬ 
ating  the  integral  is  to  fold  the  left  half  of  the  primitive  contour  in  the 
upper  half-plane  (for  x  positive)  over  onto  the  right  half  of  the  primitive 
contour,  as  shown  in  Fig.  lb  .  (A  residue  contribution  at  k  «  ik^,  where 

k„  c  2**  [m2  a"2  +  w2.  a?]*5  .  may  be  neglected  except  for  very  small  values 
1)  1  po  c  pi  iJ  *  !  1  1 

of  x  .  There  is  a  pole  at  k=0  which  gives  a  contribution  to  the  potential 
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which  is  proportional  to  sign  {x}  and  therefore  has  no  physical  signifi¬ 
cance.)  This  is  the  contour  used  by  Gould*  for  numerical  evaluation  of  the 
Fourier  inversion  integral.  It  has  the  disadvantages  that  accuracy  is  re¬ 
duced  as  x  increases,  because  of  the  rapid  oscillatory  behavior  of  the  factor 
exp(ikx),  and  that  it  docs  not  provide  for  separate  determination,  as  a 
residue,  of  the  contribution  to  the  response  arising  from  the  ion-acoustic 
pole  at  ,  which  is  dominant  uithin  a  considerable  range  of  x  . 

The  method  used  here  avoids  these  disadvantages.  It  provides  some 
additional  physical  insight  into  the  damping  of  the  wave  and  leads  to  a 
method  of  dealing  with  the  branch-cut  integrals  which  is  useful  in  the 
vastly  more  complicated  nonlinear  case. 

Hie  physical  understanding  of  Landau  damping  of  a  weakly  clamped  wave 

involves  an  interaction  between  the  wave  and  particles  with  velocities  very 

4 

near  the  phase  velocity  of  the  wave.  Presumably  the  case  of  strong  damp¬ 
ing  considered  here  involves  interactions  between  the  wave,  which  consists 
of  a  superposition  of  phase-velocity  components,  and  particles  in  the  cor¬ 
responding  band  of  phase  velocities.  Therefore  a  transformation  of  variable 
of  integration  from  k  to  the  dimensionless  complex  phase  velocity 
C  e  u>0/ka^,  which  is  considered  a  more  "natural"  variable,  is  introduced. 

The  usual  dimensionless  variables  z  =  w^x/a^  ,  f  =  w(/topi  *  *  "co^°OxO  * 

and  ¥  *  T./T  arc  introduced.  The  mass  ratio  is  p  =  m  /m.  .  'Making  use 

ci 

of  the  relation 
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the  following  separation  into  an  "ion-like"  and  an  "electron- like"  integral 
is  achieved: 


(Since  the  factor  exp(iz£~  )  has  an  essential  singularity  at  c=0  ,  giving 
z  a  small  positive  imaginary  part  assures  the  existence  of  the  integral.) 

The  dielectric  functions  are  here  and  henceforth  indicated  as  functions  of 
the  alternative  arguments  C  and  f2  ,  The  form  of  Eq.  (5)  displays 
clearly  the  role  of  ion  and  election  distribution  functions  in  determining 
the  character  of  the  response.  The  usual  factor  K"*  is  modulated  in  each 
integral  by  the  velocity  derivative  of  the  distribution  function,  which 
indicates  the  range  of  particle  velocities  which  interact  with  each  phase- 
velocity  component  of  the  wave.  The  ion-like  integral,  which  is  strongly 
affected  by  the  ion-acoustic  root  of  the  dispersion  relation  at  =  w0^ai  # 
principally  determines  the  response  from  the  grid  to  the  beginning  of  the  inter¬ 
ference  region.  See  the  dashed  curve  of  Fig.  2.  In  this  region  electrons 
easily  follow  the  low-frequency  ion  motion  and,  as  is  well  known,  neutralize 
the  ion  charge  density  quite  effectively.  In  the  region  where  the  ion-like 
and  electron-like  integrals  are  of  comparable  magnitude  the  interference  is 
obtained.  At  large  distances  the  response  is  principally  determined  by 
the  electron-like  integral;  the  phase-mixing  factor  exp(izt"*)  limits  the 
phase  velocities  involved  to  values  greater  than  those  of  most  ions  so  that 
only  an  electron  charge  density  arises.  This  part  of  the  response  corresponds 
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to  electron  shielding  of  the  disturbance  produced  by  the  grid. 

By  deforming  the  contours  so  that  they  proceed  from  the  origin  in  the 
direction  of  legative  imaginary  values,  pass  below  the  pole  at  ,  and  that 
pole  only,  and  approach  infinity  well  within  the  range  -ir/4  <  arg(c)  <  tf/4  , 
one  expresses  the  response  as  the  sun  of  the  residue  of  the  pole  at  and 

of  the  two  integrals. 

The  ion- like  integral  may  be  evaluated  along  the  path  of  steepest 

-1  2 

descents  for  the  function  exp(iz£  -  C  )  »  shown  in  Fig.  lc,  or  along  the 
simpler  contour  shown  in  Fig.  Id,  which  is  independent  of  z  .  The  path  of 
the  steepest  descents  is  the  archetype  of  desirable  contours  in  that  it 
transforms  the  dominant  exponential  behavior  of  the  integrand  into  a  gaussian; 
for  the  simpler  contour  considered,  the  exponential  becomes  predominantly 
damped  (as  opposed  to  oscillatory)  as  £  approaches  the  origin  and  infinity. 
Thus  the  oscillatory  behavior  >;iich  plagues  an  integration  along  the  real 
axis  for  large  x  is  avoided  here. 

The  evaluation  of  the  electron-like  integral  requires  a  somewhat 
different  treatment.  Since  its  integrand  decreases  slowly  as  |c|  increases, 
the  dimensionless  wave  number  is  an  appropriate  variable  and  the  contour  of 
Fig.  le  is  a  suitable  contour.  The  practical  difference  between  using  one 
variable  and  the  other  is  that  in  the  numerical  integration  one  treats  the 
variable  chosen  as  having  a  distribution  of  discrete  values  which  is  not 
radically  different  from  a  uniform  distribution.  The  path  of  steepest  descents 
for  the  function  expfizz;”1  -  p^’C2)  or  a.  contour  which  bears  the  same 
relation  to  it  as  Fig.  Id  bears  to  Fig.  lc  is  unsatisfactory  because  a  large 
number  of  roots  of  the  dispersion  relation  arc  swept  past  in  going  to  such 
a  contour  from  the  positive  real  axis. 
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Computations  were  performed  both  for  the  dipole  grid  model  and  for  the 
model  of  the  double  grid  with  finite  spacing.  The  latter  case  was  considered 
in  order  to  gain  an  appreciation  of  the  importance  of  the  assumed  theoretical 
model  of  the  grid  excitation  in  determining  the  character  of  the  response. 

Hie  factor  [sin(kx0/2)/(kx0/2) ]  introduced  into  the  integrand  of  Eq.  (1) 
gives  the  integral  for  the  case  of  spacing  between  the  grids  equal  to  x^  . 

The  potential  for  points  x  >  xQ/2  is  <Hx)  =  xQ~J  [$d(x  +  j  Xq)  -  $d(x  -  j  x^)]  , 
in  which 


(The  pole  at  k=0  gives  unimportant  contributions  to  4>(x)  as  in  the 
dipole  limit.) 

Computations  were  performed  on  the  Culler  On-Line  Computer  of  TRW 
Systems,  Inc.,  Redondo  Beach,  California.  Figure  2  shows  the  results  for 
the  dipole  limit  and  the  finite-spacing  case  (with  =  w^x^/2a^  =  4  ) 
when  T- 1  and  f 2  «  1  .  (in  the  limit  f2  «  1  one  has 
f2K±(C,f2)-7f2K±(C)  =  -C2[z;(0  +  TZ^(M1/2  T1/2  ?))  ).  Comparison  of  the 
curves  suggests  that  the  general  character  of  the  response  is  unaffected 
by  a  choice  of  reasonable  models  of  the  grid  excitation,  but  that  fine 
details  such  as  the  precise  character  of  the  interference  between  ion  and 
electron  waves  should  not  be  taken  very  seriously. 


Si  l  , 
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in.  FORMULATION  OF  TUB  NONLINEAR  PROBLEM 


Examination  of  the  (nonlinear)  Vlasov  equation  indicates  that  the  steady- 
state  response  to  time- -harmonic  grid  excitation  at  frequency  u>q  involves  . 
frequency  components  w=0,  u>^,  2 Strong  nonlinearity  involves  sub¬ 
stantial  contributions  from  a  large  number  of  harmonics  and  a  coupled  system 
which  is  probably  difficult,  if  not  impossible,  to  solve.  Accordingly  a 
solution  is  sought  for  the  case  of  weak  nonlinearity  by  introducing  a  per¬ 
turbation-series  expansion  of  the  potential  and  species  distribution  func¬ 
tions.  (This  is  equivalent  to  the  procedure  used  by  Montgomery  and  Gorman^ 
to  study  the  initial-value  problem.)  Introducing  for  convenience  the  par- 
meter  X  ,  which  is  finally  to  be  set  to  unity,  one  obtains  the  following 
expansions  to  order  X2  of  the  potential  in  the  plasma  and  the  species  dis¬ 
tribution  functions : 


(7) 


Substituting  these  expressions  into  the  Vlasov  equation  and  collecting 
like  powers  of  X  one  obtains  in  first  order  the  linearized  Vlasov  equation 
and  in  second  order  the  set 


— — ^  - _ 
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in  which  E  =  -3$/3x.  By  proceeding  to  higher  orders  one  obtains  a  hierarchy 
of  pairs  of  equations  for  higher-order  quantities.  In  each  order  the  equation 
corresponding  to  Eq.  (9)  differs  in  having  a  different  sum  of  quadratic 
conbinations  of  lower-order  terms  on  the  right-hand  side;  the  structure  of 
the  other  equation  is  the  same  as  that  of  Eq.  (10).  For  sufficiently  small 
amplitude  of  excitation  only  the  first-order  and  second-order  components  of 
the  response  need  be  considered. 

The  second-order  equations  are  Laplace-Fourier  transformed;  the  trans¬ 
form  operator  is 


The  transform-inversion  operator  is  denoted  by  T"*,  (Primed  and  double- 
primed  transform  variables,  u*  ,  w**  ,  kf  ,  and  k"  ,  are  used  below;  the  cor¬ 
responding  operators  arc  denoted  by  T' ,  T1"*,  T"  ,  and  TM  *.)  The  dis¬ 
tribution  functions  $f0(w,k,v)  arc  eliminated  from  Eqs.(9)  and  (10)  and  the 
resulting  equation  is  solved  for  6<f»(w,k),  The  electric  field  and  species 
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distribution  functions  in  linear  theory  are  expressed  as  inversions  of  the 
transforms  E(to,,,ku)  and 


•f  /to  i>'  v)  —  F(to'  V) 

TcA  )R;  J  rr*  (fcWo')U 


in  which  E^'jk’)  is  the  transform  of  the  external  field  of  the  grid  divided 
by  the  dielectric  function.  Here  w*  initially  has  a  positive  imaginary 
part.  The  velocity  integrals  which  arise  (see  Eq.  (10))  are  analytically 
continued  into  the  entire  u/k  and  u’/k*  planes.  The  deformation  of 
the  contour  of  integration  in  the  velocity  piano  around  the  pole  at  v  *u>/k 
is  the  same  as  that  in  K(w,k);  the  deformation  around  the  pole  at  v  ^lo’/k1 
is  the  same  as  that  in  K(w*  ,k*).  The  Laplace-Fourier  inversion  of  6<f,(w>k) 
is  performed..  The  result  of  these  operations  is 


S*0*>=  f— 

F  J  (v-u>KJ)  J 


in  which  E  and  f  arc  laplacc-Fouricr  transforms  of  the  field  and  dis¬ 
tribution-function  perturbations  in  linear  theory,  and  A  denotes  the  contour 

in  the  velocity  plane  deformed  as  described  above.  Making  use  of  the  no- 
2  4 

tation  U  =  q  w  /m  a  and 
a  'a  pa  a  a 


\/(r  r/)  —  ?  ^  f  i  ^  roroc</&v  ,i 
/  \ 
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in  which  Cft  =  w/ka^  and  =  u’/J^a^  one  ra3y  express  the  response  as 


Sc 


-i 


k3K(w,lO^ 


(15) 


For  a  plasma  composed  of  electrons  and  of  one  species  of  singly-charged 
r2 

ions,  Ue/U.  =  -T  ,  Performing  an  integration  by  parts  and  making  a  par¬ 
tial  fraction  expansion  in  Eq.  (14),  one  obtains 


i/fr Cz'ft'i-zV;)]  ,  z"(?) 
Vw,S)-'  (c-?os  W 


(16) 
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I V .  DOMI N ANT- PO  U;  APPR0X1 MAT1 ON 

The  practical  impossibility  of  evaluating  6$(x,t)  without  introducing 
a  further  simplification  is  now  apparent.  Application  of  the  convolution 
theorem  to  T{  )  in  Eq.  (15)  results  in  a  single  integration  with  respect  to 
k*  which  contains  a  quadratic  combination  of  a  function  of  k'  and  a  func¬ 
tion  of  (k-k1).  The  inversion  k->x  leads  to  an  integration  with  respect 
to  k  .  In  addition  to  the  requirement  of  evaluating  a  double  Fourier  inte¬ 
gral  is  .the  further  complication  that  the  functions  of  k*  and  (k-k') 
within  the  double  integral  have  branch-point  pairs  in  the  k  and  k'  planes 
when  k 1  =0  and  (l;-k,)=0,  just  as  there  is  a  branch-point  pair  at  k-0  in 


the  linear  problem.  An  escape  from  this  impasse  is  provided  by  the  dominance 
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of  the  exponentially  damped  contribution  of  tbe  residue  of  the  ion-acoustic, 
pole  over  a  wide  range  of  x  in  the  linear  problem.  The  deviation  from 
this  response  near  the  grid  is  so  strongly  damped  that  its  contribution  to 
the  nonlinear  response  is  probably  highly  localized;  the  deviation  due  to  the 
electron  wave  at  large  distances  is  at  such  a  low  level  that  it  may  be 
neglected. 

The  dominant-pole  approximation  is  obtained  by  representing  the  steady- 
state  electric  field  in  linear  theory  by  the  residue  contribution 


E(x,t )  =  -f  eA  -  P  [-  ( h 1 A I  ~  ^0t )]  +  C.c. 


(17) 


in  which 


(18) 


and  kj  =  • 

To  obtain  the  double -frequency  response,  it  is  sufficient  to  consider 
the  partial  linear  response  given  by  the  first  term  of  F.q.  (17),  The  cor¬ 
responding  nonlinear  response  plus  its  complex  conjugate  is  the  real  double- 
frequency  response.  The  Laplace-Fourier  transform  of  the  partial  excitation 
is 


!  •;  r,. 

i.:V' 


(19) 


Examining  Eqs.  (1)  and  (5)  one  notes  that  the  first  term  of  this  expression 
comes  from  K+(w0#k)  and  that  the  second  term  comes  from  K_(u>g#k).  Because 
of  the  replacement  of  the  dielectric  function  in  the  transform  of  the  response 
in  linear  theory  by  its  dominant-pole  component  it  is  necessary  to  indicate 
whether  the  plasma  dispersion  functions  in  (see  Eqs.  (15)  and  (16)) 

are  plus  or  minus  functions.  This  is  accomplished  by  the  notation 
V  which  s  and  s*  are  plus  or  minus;  the  first  subscript 

is  associated  with  the  first  argument  and  the  second  subscript  with  the 
second  argument. 

As  the  first  step  in  determining  the  explicit  form  of  Eq.  (15)  in  the 


dominant -pole  approximation  with  the  partial  excitation  described  above, 
one* obtains  the  inversion 


Here  0(x*)  is  the  unit  step  function;  equals  Cj  for  a=i 

and  for  ct~c  .  The  subscripted  parentheses  indicate  that  plasma 

dispersion  functions  of  argument  may  be  plus  or  minus  functions.  The 
further  steps  indicated  in  Eq.  (15)  are  performed;  one  obtains  6<j>  (x,t) 

=  cxp(-2ii0yt)6<}>+(x)  ,  in  which 

v  "So  '(:<)  --  iA"vi  °-0  U,( 


f'cxn  0;  U) 

Vi,.  fe,-r-3 

(i‘-R 

"  ’  ?■■  ■.)  j 

Haro  o  has  been  set  equal  to  ?u  so  that  C  -2t.^/];r'(/  ,  The  defiriilit 
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of  C|a  is  unchanged. 

TTie  folding  of  the  primitive  contour  in  the  k  plane  and  the  transforma¬ 
tion  of  the  variable  of  integration  from  k  to  £  ,  which  is  now  defined  to 
be  equal  to  2(i>0/ka^  ,  are  performed.  The  contour  in  the  £  plane  is  now 
deformed  so  that  it  proceeds  from  the  origin  in  the  direction  of  negative 
imaginar>  -  lues,  passes  below  the  pole  at  Cj  ,  and  approaches  infinity 
well  within  the  range  -n/4<arg{c}<if/4  .  Such  a  contour  is  denoted  by  C  . 

As  in  the  linear  problem  there  is  a  simple  pole  at  k=0  ,  which  has  no  physi¬ 
cal  significance.  For  the  case  f*«l  (f  =  “(/“pi*  as  before)*  t0  which  con“ 
sideration  is  now  limited,  f2K+(2w0,k)  -  j.  f2K±(S)  .  For  purpose  of  com¬ 
parison  with  integrals  in  the  linearized  problem  it  is  convenient  to  consider 
[6*+(x)/£6  (-2itiC.)],  in  which  C.  =  Cja=i  and 


9  JIIWIVl 


^  £  ■'>.rr,iL. 
f'  r~  %  rrr  *  -  ■• 
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in  which 


X 


yd)  __  f  _ JtexpCS-gS"*) _ 

1  -  Jc  (-e*te-sBeK*meK<® 

|fr  K.fe!  \4+  fe/i)  -  fa  K+  K)V-f  (%)  ] 
(C-Q 

ff?“  K-C0Vf-(^)--FgK4(c)\/--^)] 
(c+&) 


I 


(e) 


X 


f  d?e'4p( 

A  1  •p'm/i  \ 

I  c-eiri)  r*-~| 

/N  * 

]CfsK-(?a 

fr-Sk) 


:sa>  -f  Ei 
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- 
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(24) 


(25) 


and  R  is  the  residue  contribution,  which  will  be  discussed  and  evaluated 

1/21/2  1/2  1/2 

below.  For  compactness  the  notation  £  =  y  '  T  C  and  Cj  =  V  VT  '  Cj 

has  been  introduced.  The  definition  of  z  remains  unchanged,  z  =  t^x/a^  . 

Expression  of  1^  and  1^  as  sums  of  ion- like  and  electron-like 

integrals,  that  is  integrals  whose  integrands  have  the  exponential  behavior 
-12  -1-2 

cxp(izC  -  C  )  and  exp(izC  -  pT£  )  ,  respectively,  is  achieved  by 
manipulations  which  are  here  described  briefly.  By  expressing  ZM CC)  *n 
terms  of  Z'(0  in  Eq.  (16),  one  obtains 


V(C,C')={  (C-CT*  Z’M  +  £| ?(?-«]"} 


I  imWHMI  iiwiiaifWWIilWiMli^^  i 


In  1^  and  1^  only  the  dependence  of  the  V's  remains.  To 


indicate  the  structure  of  the  manipulations  one  may  denote  the  V*s  by 


W  s  -  W  s  V-t^-ho)  •  in  whi 


which 


iAau= p&i+Q(c<)z±fc). 


The  P's  and  Q's  depend  on  the  species  a  and  on  the  suppressed  ±  index 


associated  with  the  arguments  ±  5^  .  But  they  do  not  contain  Z'+(^ 


and  hence  are  independent  of  the  ±  index  associated  with  ^  .  Appropriate 


re) 

elements  of  the  integrands  of  Iv  *  and  Iv  9  are  expressed  as 


[f  eK-(r)  w+(rj  -  f 1 sK+(?)  iv-(Q] 


= V4&)  [f e  kGM'K*®] +Q  ( d  [z.fe!  -z/c  <y] . 


The  desired  expression  of  1^  and.  1^  as  sums  of  ion-like  and 


electron- like  integrals  is  achieved  by  making  use  of  liq.  (4). 


Defining,  for  n  =  2,3,  the  functions 


and  making  use  of  the  relation  Z’4(- C)  s  Z'^(C)#  one  obtains  finally 


for  a  e  i,c 


-i  w  K// 


ts,  -Ji. 


in  which 


(e<f3M_  s  f  e*  P  (2i  h  r*)G (a;3n)  Jr 


ir  AJ  c-e  p«K^)if ^)] 
c 


The  func  'ns  G^n^  are 


.ftlS) 


e*p  wr>[Zt®  +t  z+Ysi]  d3(  r) 


G8is)=exp(-fe)c*[e-T(i-efe)ZY©]D£(y 


rC£es)  4-^s.  . ,  .  ,| 

G  =H*T  expf-?5)?"' 


G{l!“LFZ  j'expGr-^  fe[2  +  (1-3  r£)z+'(s]pe 


Qfeli~— pTJf  eypt-^^fz^-Z/O]  D,(g 


In  Eqs.  (36)  and  (38)  the  relation  D3(£)  =  y'3/2T~3/2D3(c)  has  been 

•  1  AS  1 

used;  in  Eqs.  (37)  and  (39)  the  relation  ^(C)  =  y  T  02(C)  has  been 
used. 
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V.  SIMPLI FIGVTION  AND  NUMERICAL  RESULTS 


The  numerical  evaluation  of  the  branch-cut  integrals  nay  be  simplified 

considerably  by  the  introduction  of  approximations  and  an  ordering  of  terns 

based  on  the  smallness  of  1 i1^2  .  In  lowest  order  1^  may  be  neglected. 

In  next  order  cne  element  of  1^  ,  namely  I^**2^  ,  makes  a  contribution 

to  the  result  which,  although  small,  is  enhanced  by  the  character  of  its 

integrand  above  that  indicated  by  the  ordering. 

The  relation  [Z * ^ +  fz*  +  (^)]  =  0  implied  by  the  dispersion 

relation  facilitates  the  comparison  of  the  relative  magnitudes  of  the 

corresponding  elements  of  1^  and  3^  which  contain  .  Using 

this  relation  one  obtains  -  p  and  [G^e^/G^ee^]  =  pT”1 

The  corresponding  elements  of  and  1^  which  contain  D2  are 

(ie2)  fee21 

now  compared.  Consider  first  Iv  and  lv  .  In  the  small -argument 

region  of  the  contour  C  ,  |tj  ~  4  ,  the  factor  [2  ♦  (1-2C  )Z'+(C)]  , 

which  appears  in  G^ie2^  ,  and  the  factor  [2T  -  (1-2£2)Z'  +  (C)]  ,  which 
(ee2) 

appears  in  Gv  ,  are  both  of  order  urity.  Thus  the  ratio  of  the 

absolute  value  of  the  small -argument  part  of  I^e2^  to  that  of  I^ee2^ 

is  approximately  p  .  In  the  asymptotic  region  |c|  £  4  the  two  factors 

have  the  behavior  [2  +  (l-2?2)Z'+(0]  [=  CZ’  +  (  )]  'v  -  2c"2  ami 

[2T  -  (l-2£2)Z'+(s) ]  2T  .  Thus  the  ratio  of  the  absolute  value  of  the 

f ie2)  fee  21 

asymptotic  part  of  I  to  that  of  Iv  is  less  than  p  . 

The  last  comparison  to  be  made  is  that  of  I^*2^  and  I  ^e*2^  .  in 
the  range  |c|  ~  4  over  which  numerical  evaluation  of  ion- like  integrals 

r 1  ? 

is  performed,  the  corresponding  factors  [2  -  T(J-2£  )Z'+(£)]  ,  which 
appears  in  G^*2^  ,  .and  CZ"+(£)  ,  which  appears  in  G^c*2^ 


,  are  both  of 
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order  unity.  Therefore  the  ratio  of  the  absolute  value  of  to 

that  of  I^c*2^  is  0(p*^2)  .  It  will  be  seen  that,  due  to  the  absence 
of  inverse  powers  of  (C-^)  in  the  lowest  order  of  an  expansion  of  the 
integrand  of  the  sura  +  j(e*2)]  ,  I^**2^  contributes  more 

importantly  to  the  total  response  than  the  ordering  indicates.  Therefore 
it  is  retained. 

The  integrands  which  have  been  retained  are  now  simplified.  Consider 
first  the  sum  of  ion-iike  integrals,  [I^ei3^  +  I^e*2^  +  i^ii2)]  #  Expanding 
the  factor  [Z*  +  (£p  -  Z'+(£)]  which  appears  in  .  G^ei3^  in  powers  of 
(S-tj)  and  approximating  Zu+(£)  in  G^ei2^  by  ,  one  obtains 

in  lowest  order 


[Gsw+G^»]=1r5f'e,p(-f^3[£f1z+'fe)](f-rirf 


The  reason  for  retaining  i^*2^  should  now  be  apparent.  Although  G^ii2^ 

is  formally  smaller  than  [G^ei^  ♦  G^ei2^]  by  the  factor  p1^2  ,  the 

-2 

presence  of  a  factor  (C-C^)  in  the  former  and  the  absence  of  negative 
powers  of  (C-C^)  in  the  latter  result  in  a  contribution  of  i^2^ 
substantially  larger  than  that  indicated  by  the  ordering.  (The  next  order 
of  the  expansion  of  ♦  G^e*2^]  makes  a  negligible  c  .  trjbution  to 

the  result,  despite  the  fact  that  one  of  two  terms  contains  a  factor 
(C-Cj)  .)  Since  only  lowest  order  terms  in  p  '  are  retained  (with  the 
exception  of  the  anomalously  large  contribution  I^2^  )  the  approximations 
Z’^CSj)  et  -2u1/2i  and  Z"+(£)  «  -2  may  be  made  in  Eqs.  (33)  and  (40), 
respectively. 

fce3)  (002") 

Consider  now  the  sum  of  electron- like  integrals,  [Iv  ♦  Iv  ]  . 


Tfcey  are  evaluated  as  the  sun  of  two  parts:  a  snail-argent  component 
evaluated  along  the  sane  contour  as  the  ion-like  integrals  and  a  large- 
argument  conponent  evaluated  along  a  contour  in  the  complex  phase  velocity 
plane  of  the  same  general  character  as  that  of  Fig.  le  but  with  an  upper 
limit  which  is  the  reciprocal  of  the  upper  limit  of  integration  in  the 
C  plane.  In  the  small-argument  region  an  expansion  in  p1/2  gives 


We  term  of  order  Vl,?  makes  an  anomalously  large  contribution  to  the 

integral  because  [Z-+(«  .  ft  v.xnlshes  at  a  point  with  a  distance  of 

order  p  from  ^  and  is  therefore  retained.  I„  the  large-argument 

region  Z-+(C)  is  represented  by  a  suitable  number  of  terms  of  its 
asymptotic  expansion. 

The  residue  contribution  to  Hq.  (23)  is 


Wirriimiiyiinri  unmn  "i;rm>myr!i';ni.'ii  irpwiSHfiP 


One  expands  the  factors  of  Eq.  (42)  in  powers  of  (C~Cp  to  obtain  an 

explicit  expression  for  the  residue,  .  To  develop  an  approximate  expression 

which  is  convenient  for  computational  purposes  one  expands  the  residue  in 
1/2 

powers  of  p'  .  The  final  result  is 


‘L  *r» 


R  =  |IS  T  U2K+fe  2"  expfcizrCso  +  SiM 


in  which  sQ  =  -  (ir1/2/2)i  +  0(p1/2)  and  s,  =  p1/2  T_3/2(-i/2)  [Z'^Cq) 

\/2  "2 

♦  2tt /z  r  i]  .  The  term  proportional  to  z  is  retained  even  though  it 

1/2 

is  formally  of  order  V  because  it  can  make  a  substantial  contribution 
to  the  result  at  large  values  of  z  which  are  within  the  probable  range 
of  validity  of  the  dominant -pole  approximation  for  1  T£t  values  of  the 

ratio  T  /T.  . 

©  1 

The  results  of  numerical  calculations  are  shown  in  Tigs.  3  and  4  for  the 
cases  T  *=  1,  .5,  and  .25  .  It  is  difficult  to  estimate  the 

range  of  z  over  which  the  dominant-pole  approximation  is  valid.  If  the 
dominant-pole  component  of  the  linear  response,  Eq.  (17),  is  cciifidered 
to  be  mixed  with  itself  in  the  absence  of  further  dispersion  or  damping, 
the  double- frequency  harmonic  is  attenuated  as  expf-21ro(}^;xj  .  The  change 
in  x  dependence  x  *  2x  is  adopted  as  the  basis  for  choosing  a  maximum 
value  of  z.  For  the  equal  temperature  case  the  deviation  of  the  linear 
response  from  the  dominant -pole  approximation  is  small  (except  very  near 
the  grid)  until  z  -  22  .  Thus,  z-  11  is  chosen  us  the  maximum  value 
for  the  calculation  of  the  double- frequency . response.  For  the  case  Te/]\  -  2, 
the  corresponding  maximum  value  for  the  linear  response  is  z  -  38.  (See 
the  results  of  Gould,  Fig.  G  of  Ref.  1.)  Accordingly,  z  =  19  is  selected 
as  the  maximum  value  for  the  calculation  of  the  double- frequency  response. 


****** 


For  this  and  higher  values  of  the  temperature  ratio  the  double- frequency 
response  has  a  characteristic  repetitive  pattern;  accordingly,  in  the 
interest  of  conserving  computer  time  the  corresponding  increases  in  the 

maxima  values  of  z  are  not  made  beyond  T  /T.  *=  2  . 

6  1 

The  principal  features  of  the  numerical  results  are  three:  a 

repetitive  modulation  of  the  spatially  damped  behavior  of  the  logarithm 

of  the  absolute  value  of  the  potential  perturbation;  a  reduction  in  spatial 

damping  as  the  electron-to-ion  temperature  ratio  increases;  and  a  transition 

in  the  overall  spatial  rate  of  change  of  phase  with  temperature  ratio 

increase,  from  a  small  rate  characteristic  of  a  non-wave- like  disturbance 

similar  to  the  Electron  wave*'  at  large  distances  in  the  linear  problem, 

to  a  large  linear  rate  characteristic  of  the  mixing  of  the  dominant-pole 

coriponent  with  itself  and  similar  to  the  behavior  of  the  phase  of  the  ion- 

acoustic  wave  in  the  linear  problem.  The  modulation  of  the  amplitude  of 

the  response  and  the  non-wave-like  character  of  the  response  for  values  of 

the  electron-to-ion  temperature  ratio  close  to  unity  are  consequences  of 

the  fact  that  the  electron- like  integral  in  the  double- frequency  response 

is  not  diminished  relative  to  the  ion-like  integral  by  the  presence  of  a 
1/2 

factor  p  ,  as  is  the  case  in  the  linear  response. 

VI.  ZERO -FREQUENCY  RESPONSE 

The  nonlinear  response  at  zero  frequency  is  now  shown  be  a 

polarization  of  the  plasma  with  no  associated  current  den  r. .  This 

result  is  established  in  the  context  of  the  assumptions  e  1 1-amplitude  . 
excitation,  uniform  and  stationary  distribution  function  *  .bscnce 
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of  excitation,  and  one-dimensionality.  There  are  a  number  of  conditions 
in  the  experiments  which  are  inconsistent  with  these  assumptions;  the 
inclusion  of  these  complcations  would  make  the  analysis  vastly  more  difficult. 

The  zero- frequency  component  of  potential  in  the  plasma  must  be  an  even 
function  of  x  because  there  is  no  preferred  direction  in  space  when  an 
average  over  one-  cycle  of  the  excitation  is  considered.  The  only  possible- 
source  of  a  zero- frequency  component  of  potential  in  this  theory  is  the 
inequality  of  the  charge-to-mass  ratio  for  the  two  species,  which  would 
give  a  potential  which  is  an  even  function  of  x  . 

The  Vlasov  equation  implies  species  continuity  equations 


ltn«  +  li1  r«=o  <«> 

in  which  is  the  particle  current  density  of  species  a  .  For  the 
zero- frequency  component  this  equation  states  that  the  species  particle 
current  densities  are  divergence less.  The  spatial  symmetry  of  the  zero- 
frequency  component  of  the  potential  implies  that  species  particle  current 
densities  are  zero  at  x  =  0  ;  hence  they  must  be  zero  everywhere. 

The  zero- frequency  component  of  the  lowest  order  nonlinear  contribution 
to  the  potential  was  determined  in  Ref.  6.  Since  the  zero- frequency 
component  is  probably  difficult  to  detect  in  an  experiment,  that  development 
is  not  included  here. 


VII.  CONCLUSIONS 


The  double- frequency  response  produced  by  the  nonlinear  interaction  of 

a  grid-excited  ion-acoustic  wave  with  itself  exhibits  complicated  behavior 

for  values  of  the  electron-to-ion  temperature  ratio  near  unity.  The  electron- 

wave  component  of  the  response  is  not  diminished  by  the  presence  of  a  factor 
1/2 

p  ,  as  occurs  in  the  linear  case.  This  indicates  the  operation  of  an 
electron  shielding  mechanism  in  the  double  frequency  response.  As  the 
temperature  ratio  approaches  a  value  of  four,  the  response  comes  closer  to 
displaying  exponentially  damped  behavior. 

The  complications  involved  in  determining  nonlinear  response  to  grid 
excitation  are  reduced  considerably  by  expressing  the  Fourier  inversion 
integrals  which  result  as  sums  of  "ion-like"  and  "electron-like"  integrals. 
Physical  interpretation  of  the  results  of  a  nonlinear  grid  excitation  problem, 
is  rendered  somewhat  less  difficult  by  such  a  representation. 
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Figure  Captions 

Fig.  1.  Integration  contours  for  Fourier  inversion  integrals. 

(a)  Primitive  Fourier  inversion  contour  in  k  plane  and  mapping 

of  branch- cut  and  branch  points  of  Z  +(Cj  for  0)  =  +  ie  . 

(b)  Folding  of  left  half  of  primitive  contour  onto  right  half  of 

primitive  contour  in  k  plane.  First  few  members  of  infinite  set 

of  roots  of  dispersion  relation  K+  =  0  ai'e  indicated. 

-1  2 

(c)  Path  of  steepest  descents  for  function  exp(iz£  -  C  )  and 
mapping  of  roots  of  dispersion  relation  shown  in  (b)  in  £  plane. 

(d)  Simple  contour  for  evaluation  of  ion-like  integrals. 

(e)  Contour  in  k  plane  for  evaluation  of  electron-like  integrals 
which  reduces  oscillatory  behavior  of  integrand. 

Fig.  2.  Numerical  results  in  linear  theory  (cesium  plasma). 

Fig.  3.  Numerical  results  in  nonlinear  theory  (cesium  plasma):  natural 

logarithm  of  absolute  value  of  [I^‘  +  I  ^  +  R]  as  a  function  of 
z  for  temperature  ratios  T  =  1,  .5,  .25  . 

Fig.  4.  Numerical  results  in  nonlinear  theory  (cesium  plasma):  argument 

of  [1^  +  1^  +  R]  as  a  function  of  z  for  temperature  ratios 


T  *  1,  .5,  .25  . 
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